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Abstract — This paper presents a new approach to solving the 
short-term unit commitment problem using an improved 
Particle Swarm Optimization (IPSO). The objective of this 
paper is to find the generation scheduling such that the total 
operating cost can be minimized, when subjected to a variety 
of constraints. This also means that it is desirable to find the 
optimal generating unit commitment in the power system for 
the next H hours. PSO, which happens to be a Global 
Optimization technique for solving Unit Commitment 
Problem, operates on a system, which is designed to encode 
each unit's operating schedule with regard to its minimum 
up/down time. In this, the unit commitment schedule is coded 
as a string of symbols. An initial population of parent 
solutions is generated at random. Here, each schedule is 
formed by committing all the units according to their initial 
status ("flat start"). Here the parents are obtained from a pre- 
defined set of solution's i.e. each and every solution is adjusted 
to meet the requirements. Then, a random decommitment is 
carried out with respect to the unit's minimum down times. A 
thermal Power System in India demonstrates the effectiveness 
of the proposed approach; extensive studies have also been 
performed for different power systems consist of 10, 26, 34 
generating units. Numerical results are shown comparing the 
cost solutions and computation time obtained by using the 
IPSO and other conventional methods like Dynamic 
Programming (DP), Legrangian Relaxation (LR) in reaching 
proper unit commitment. 

Index Terms — Unit Commitment, Particle Swarm 
Optimization, Legrangian Relaxation, Dynamic Programming 

I. Introduction 

In power stations, the investment is pretty costlier and 
the resources in operating them are considerably becoming 
sparse of which focus turns on to optimizing the operating 
cost of the power station [1]. The demand for the electricity 
is varying in a daily and weekly cyclic manner and this 
creates a problem for the power system in deciding the best 
way to meet those varying demands [2]. The demand 
knowledge in the future is the main problem of the 
planning. With an accurate forecast, the basic operating 
functions like thermal and hydro thermal UC, economic 
dispatch, fuel scheduling and unit maintenance can be 
performed efficiently [3]. The energy production cost 
varies considerably between all the energy sources 
available on a power system. Moreover it requires a tool in 
order to balance the demand and generation and also to 
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dispatch the generation optimally in the most economical 
way [4]. 

Power system operators have to face the decision-making 
problems extensively because of these difficulties. The 
scheduling of generators in a power system at any given 
time is one of the decision making problems. Running all 
the units is not economical for a power system which is 
required to satisfy the peak load during low load periods 
[5], The main objective is to reduce the power generation 
costs when meeting the hourly forecasted power demands. 
UCP is the method of finding an optimal turn on and turn 
off schedule for a group of power generation units for each 
time window over a given time horizon [10]. The UCP is a 
vital area of research which concerns more attention from 
the scientific community because of the fact that even 
small savings in the operation costs for each hour can lead 
to the major overall economic savings [6]. In order to 
decide which of the available power plants should be 
involved to supply the electricity, the best choice is the UC 
which is called as Unit Commitment [7]. UCP is an area of 
production scheduling which is related to decide the 
ON/OFF status of the generating units during each interval 
of the scheduling period. This is done in order to meet the 
system load and reserve requirements and minimum cost 
which are exposed to many types of equipment, system and 
environmental constraints [8]. The UC should 
simultaneously reduce the cost of the system production 
when it satisfies the load demand, spinning reverse, ramp 
constraints and the operational constraints of the individual 
unit [9], 

Some of the recent research works related to solving the 
UCP is discussed as follows. A novel approach for solving 
the short- term commitment problem using the genetic 
algorithm based tabu search method with cooling and 
banking constraints was proposed [2]. The main aim of his 
work is to find the generation scheduling so that the total 
operating cost can be reduced, when subjected to a variety 
of constraints. A two- layer approach for solving the UCP 
in which the first layer utilizes a Genetic Algorithm to 
decide the on\off status of the units and second layer 
utilizes an Improved Lambda- Iteration technique to solve 
the Economic Dispatch Problem was presented [11]. Their 
approach satisfies all the plant and system constraints. UCP 
for four-unit Tuncbilek thermal plant which was in 
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Kutahya region in Turkey, was solved [3] for an optimum 
schedule of generating units based on the load data 
forecasted by using the conventional ANN (ANN) and an 
improved ANN model with Weighted Frequency Bin 
Blocks (WFBB). They have used Fuzzy Logic (FL) 
method for solving the UCP. Three evolutionary 
computation techniques, namely Steady-State Genetic 
Algorithms, Evolutionary Strategies and Differential 
Evolution for the UCP have compared [12]. Their 
comparison was based on a set of experiments conducted 
on benchmark datasets as well as on real-world data 
obtained from the Turkish Interconnected Power System. 
An algorithm which is used to solve security constrained 
UCP with both operational and power flow constraints 
(PFC) have been proposed [9] for planning a secure and 
economical hourly generation schedule. Their algorithm 
introduces an efficient unit commitment (UC) approach 
with PFC which obtains the minimum system operating 
cost which satisfies both unit and network constraints when 
contingencies are included. 

Alternative strategies with the advantages of Genetic 
Algorithm for solving the Thermal UCP and in addition to 
these they have developed [8]. Parallel Structure to handle 
the infeasibility problem in a structured and improved 
Genetic Algorithm (GA) which provides an effective 
search and therefore greater economy. A novel approach 
for solving the Multi- Area Commitment problem using an 
evolutionary programming technique has proposed [13]. 
Their technique was used to improve the speed and 
reliability of the optimal search process. A dynamic 
programming algorithm for solving the single-UCP (1- 
UCP) which efficiently solves the single-unit economic 
dispatch (ED) problem with ramping constraints and 
arbitrary convex cost functions have proposed [14]. The 
analytical and computational necessary & sufficient 
conditions to determine the feasible unit commitment states 
with grid security constraints have presented [15]. The 
optimal scheduling of hydropower plants in a hydrothermal 
connected system has considered [16]. In their model they 
have related the amount of generated hydropower to 
nonlinear traffic levels and also have taken into account the 
hydraulic losses, turbine- generator efficiencies as well as 
multiple 0-1 states associated with forbidden operation 
zones. 

From the surveyed research works it can be understood 
that solving the UCP gains high significance in the domain 
of power systems. Solving the UCP by a single 
optimization algorithm is ineffective and time consuming. 
Hence, we are proposing a UCP solving approach based on 
improved PSO which provides an effective scheduling with 
minimum cost. The proposed approach solves the UCP 
with less time consumption rather than the approaches 
solely based on a single optimization algorithm. 

II. Particle Swarm Optimization 

PSO first introduced by Eberhart and Kennedy in 1995 
which uses the natural animal's behavior like bird flocking, 



fish schooling, and swarm theory to produce the best of the 
characters among comprehensive old populations [17-20]. 
The basic PSO algorithm can be described as follows: Each 
particle in the swarm represents a possible solution to the 
existing optimization problem. During PSO iteration, every 
particle accelerates independently in the direction of its 
own personal best solution which is found so far, as well as 
the direction of the global best solution which is discovered 
so far by the other particles. Therefore, if a particle finds a 
promising new solution, all other particles will move closer 
to it, exploring the solution space more thoroughly [22]. 
Computation in PSO is based on a population (swarm) of 
processing elements called particles in which each particle 
represent a candidate solution [18]. The PSO algorithm 
depends on the social interaction between independent 
particles, during their search for the optimum solution [19]. 
PSOs are initialized with a population of random 
solutions and search for optima by updating generations 
[21]. All particles have fitness values which are estimated 
by the fitness function to be optimized, and have velocities 
which direct the flying of particles [20]. The fitness 
function is evaluated for each particle in the swarm and is 
compared to the fitness of the best previous result for that 
particle and to the fitness of the best particle among all 
particles in the swarm [24]. The algorithm iterates by 
updating the velocities and positions of the particles, until 
the stopping criteria is met [25]. The position (i.e. solution) 
of every individual particle will be attracted stochastically 
towards their related best positions (i.e. best solutions) in 
multidimensional solution space [20]. The PSO algorithm 
is becoming very popular due to its simplicity of 
implementation and ability to quickly converge to a 
reasonably good solution. Nowadays, PSO algorithm is 
effectively applied in power system optimization, traffic 
planning, engineering design and optimization, and 
computer system etc [23]. 

III. Problem Formulation 

The main aim is to find the generation scheduling so that 
the total operating cost can be reduced when it is exposed 
to a variety of constraints [26]. The overall objective 
function of the UCP is given below, 

T N R 

frH (Fa ( p it Pu + s it v it ) — (i) 

t~\ M h 

Where 

Ut ~ unit i status at hour t=l(if unit is ON)=0(if unit is 

OFF) 

V; t ~ unit i start up / shut down status at hour t =1 if the unit 

is started at hour t and otherwise. 

F T ~ total operating cost over the schedule horizon (Rs/Hr) 

Sh ~ start up cost of unit i at hour t (Rs) 

For thermal and nuclear units, the most important 
component of the total operating cost is the power 
production cost of the committed units. The quadratic form 
for this is given as 
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Fit(Pit)=A l P 2 l t+B i P lt+ C i 



Rs 



(2) 



Ai, Bj, Q ~ the cost function parameters of unit i (Rs./MWiir, Rs./MWhr, 
Rs/hr) 

F it(P it ) - production cost of unit i at a time t (Rs/hr) 

P i t ~ output power from unit i at time t (MW) 

The startup value depends upon the downtime of the unit. 
When the unit i is started from the cold state then the 
downtime of the unit can vary from a maximum value. If 
the unit i have been turned off recently, then the downtime 
of the unit varies to a much smaller value. During the 
downtime periods, the startup cost calculation depends 
upon the treatment method for the thermal unit. The startup 
cost Sjt is a function of the downtime of unit i and it is 
given as 



S it = So i 



1 - D t exp 



Tdown 



V 



I J 



E t Rs 



(3) 



Where 

So; - unit i cold start — up cost (Rs) 

Dj, Ej ~ start — up cost coefficients for unit i 

A. Constraints 

Depending on the nature of the power system under 
study, the UCP is subject to many constraints, the main 
being the load balance constraints and the spinning reserve 
constraints. The other constraints include the thermal 
constraints, fuel constraints, security constraints etc. [26] 

B. Load Balance Constraints 

The real power generated must be sufficient enough to 
meet the load demand and must satisfy the following 
factors given in (4). 



ffiPn 



pa 



(4) 



Where 

PD , ~ system peak demand at hour t (MW) 

N ~ number of available generating units 

U(0,1) ~ the uniform distribution with parameters and 1 

UD(a,b) ~ the discrete uniform distribution with parameters 

a andb 

C. Spinning Reserve Constraints 

The spinning reserve is the total amount of real power 
generation available from all synchronized units minus the 
present load plus the losses. It must be sufficient enough to 
meet the loss of the most heavily loaded unit in the system. 
This has to satisfy the equation given in (5). 

2>max,t/,., >=(PD t +R,);\<t<T (5) 

;=i 

Where 
Pmaxj ~ Maximum generation limit of unit i 



R t ~ spinning reserve at time t (MW) 
T ~ scheduled time horizon (24 hrs.) 

D. Thermal Constraints 

The temperature and pressure of the thermal units vary 
very gradually and the units must be synchronized before 
they are brought online. A time period of even 1 hour is 
considered as the minimum down time of the units. There 
are certain factors, which govern the thermal constraints, 
like minimum up time, minimum down time and crew 
constraints. 

Minimum up time: 

If the units have already been shut down, there will be a 
minimum time before they can be restarted and the 
constraint is given in (6). 

Tor\ >Tup t (6) 

Where 

Tonj ~ duration for which unit i is continuously ON (Hr) 

Tup i ~ unit i minimum up time (Hr) 

Minimum down time: 

If all the units are running already, they cannot be shut 
down simultaneously and the constraint is given in (7). 

Toff>Tdo Wt} (7) 

Where 

T down ; ~ unit i minimum down time (Hr) 

T off ; ~ duration for which unit i is continuously OFF (Hr) 

E. Must Run Units 

Generally in a power system, some of the units are given 
a must run status in order to provide voltage support for the 
network. 

F. Ramping Constraints 

If the ramping constraints are included, the quality of the 
solution will be improved but the inclusion of ramp-rate 
limits can significantly enlarge the state space of 
production simulation and thus increase its computational 
requirements. And it results in significantly more states to 
be evolved and more strategies to be saved. Hence the CPU 
time will be increased. 

When ramp-rate limits are ignored, the number of 
generators consecutive online/offline hours at hour t, 
provides adequate state description for making its 
commitment decision at hour (t+1). When ramp-rate limits 
are modeled, the state description becomes inadequate. An 
additional status, generators energy generation capacity at 
hour t is also required for making its commitment decision 
at hour (t+1). These additional descriptions add one more 
dimension to the state space, and thus significantly increase 
the computational requirements. Therefore, we have not 
included in this algorithm. 

IV. Improved PSO Algorithm for solving UCP 

The proposed IPSO Algorithm is to determine the units 
and their generation schedule for a particular demand with 
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minimum cost. In this manner, with the assistance of PSO 
we determine the same for different possible demands. The 
problem is divided into two stages; one is for determining 
the minimum cost for a particular demand and another is 
for determining the minimum cost for unit commitment 
during all the periods. But the demand varies during all the 
periods. Hence, different possible demands are need to be 
generated which can be performed by the IPSO algorithm. 
PSO is used to determine the optimal generation schedule 
for a particular demand. The steps of the algorithms which 
is used for our approach is demonstrated in the "Fig. 1". 

As depicted in "Fig. 1", for a power demand of P^ , 

initially, a population of random individuals is taken. The 
random individuals include random particles and their 
velocities. Hereby, a logical algorithm is utilized to 
generate the initial random solutions of particles which can 
be discussed as follows 

1 . Generate an arbitrary integer r which satisfies the 
condition r < n . 

th 

2. For the r unit, generate a random integer 

indicating the power generated by the unit which 
should essentially satisfy the condition 

n(min) ^ r> < n(max) 

Sr ~ g, ~ Sr 

3. The remaining power to be generated i.e. 
P d — P is subjected for the following decision, 



P d = 



g. 



g. 



if P d ~Pg<P„ 



th 



else 



(8) 



4. Allot P d to the next unit to generate (now, let the 

next unit as r ) whose maximum limit of power 
generation is greater than the remaining free units. 

The allotment of Pa to the unit is based on the 



following condition 
■Jf 

■Jf 

;tf 



p g ,= 



3 (min) 



n(min) < p < D 
g, ~ d — g r 

(min) 



(max) 



p(max) 

Sr 



P^Pg, 

P -^ D( max ) 

r d > r g r 



(9) 



5. Determine P d -P . lfP d -P <0, go to step 

1; IfP d -P g > 0, go to step 3; lfP d - P g = 0, 

then terminate the criteria. 
By the above mentioned algorithm, a vector is obtained 
which represents the amount of power to be generated by 
each unit. Hence, some different possible vectors are 
generated by repeating the algorithm and it can be given as 



P = 



3(0 



p(i) 
82' 



>(0 



(10) 



Equation (13) represents the initial particles that are 
satisfying the constraints given in equation (4) and (5) are 
generated as random initial solutions for the PSO 
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algorithm. In parallel, random velocities are also generated 
for the corresponding particles as follows 



(0 (0 

'] > V 2 > 






(11) 

In equation (14), the velocities for each particle element are 
randomly generated within the maximum and minimum 

limit and so the each element of the velocity vector v ; - 

satisfies v mm < v- ( j) < v max . After determining the 

initial particles and their corresponding velocities, the 
particles are evaluated by the evaluation function which is 
given by 

min£(/^,)^ + S^,A 
i=i h 



(12) 



Generate random particles P^ and 
elo cities V: 




^Determine Pi . £s^ s. 

best c -&es 




Gtitain the optirtvuuri 

generation schedule 



(13) 



Figure 1. Steps involved in PSO to determine the optimal generation 
schedule 

Based on the evaluation function given in equation (15), 
Pfo est and g b est for the initial particles are determined. 
Then new velocities are determined as 
vr(j) = w * vr > U) + Ci * ai * 

[p bes( u)-pr'(m+c 2 *a 2 *[g best u)-prv)] 

where, 1 < i < I , l< j < n , V; (j) stands for current 

velocity of the particle, V ( - (j) stands for new velocity 
of a particular parameter of a particle, fljand a.j_ are 
arbitrary numbers in the interval [0,1] , Cj and C 2 are 
acceleration constants (often chosen as 2.0) andw is the 
inertia weight that is given as 
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li' 



w = w„ 



max 



■ w. 



min 



x/ 



(14) 



^max . 

^max 

where, W,_ ,,„ and W~.,„ are the maximum and minimum 

7 mdx iiiiii 

inertia weight factors respectively which are chosen 
randomly in the interval [0,1] , / is the current number of 

iteration and / max is the maximum number of iterations. 

The velocity of such newly attained particle should be 
within the limits. Before proceeding further, this would be 
checked and corrected. 

v „ e w ( . ) = [ V maxO') J V? eW (j)>V rmx (j) 

UmnO") ;v? ew U)<v min U) 

Depend upon the newly obtained velocity vector, the 

particles are updated and obtained as new particles as 

follows 

P^ U) = p new (j) + v new U) (lfi) 

Then the parameter of each particle is checked whether it 
is ahead the lower and upper bound limits. The minimum 
and maximum generation limit of each unit is referred by 
the lower and upper bound values respectively. If the new 
particle infringes the minimum and maximum generation 
limit, then a decision making process is performed as 
follows 



\j)- 



p(max) 

8i 
p(min) 



I if 
I if 



gi 



v ( i) > p( max ) 

(min) 



(17) 



P new ( j) < P 

Si KJJ Si 



The newly obtained particles are evaluated as mentioned 
earlier and so P\, est for the new particles are determined. 

With the concern of P\, est and thcgf, est , newgj, es( is 

determined. Again by generating new particles, the same 
process is repeated until the process reaches the maximum 
iteration / max . Once the iteration reaches the / max , the 

process is terminated and so that a generation schedule of 
all the units with minimum cost is obtained which will 
meet the demand at the particular period. In the similar 
fashion, the optimum generating schedule for all the 
possible demand set is determined. So, a complete training 
set which includes the various possible demands and the 
corresponding optimum generation schedule is generated. 

V. Results and Discussion 

The proposed intelligence technique for UCP which is 
based on the IPSO has been implemented in the working 
platform of MATLAB (version 7.8). We have considered 
an Indian thermal power system with seven unit's utility 
system for a time span of 24 hours for evaluating the 
performance of the proposed technique. The operation data 
for the system is given in the Table I. The daily load data of 
10, 26, 34 unit systems are shown in Table II. 
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Table I 
Operation Data Of Seven Units Utility System 



Unit 




Running Cost 




Start-up 


ost 


Huh 

(MW) 


(MW) 


Ai Bi d 
(RsMWlf) (RsMWh) (Rs) 


So; 
(Rs) 


Di 
(Rs) 


Ei 
(Rs) 


1 


0.255 


70 


750 


4250 


29.5 


10 


15 


60 


2 


0.198 


75 


1250 


5050 


29.5 


10 


20 


80 


3 


0.198 


70 


2000 


5700 


28.5 


10 


30 


100 


4 


0.191 


70 


1600 


4700 


32 5 


9 


25 


120 


5 


0.106 


75 


1450 


5650 


32 


9 


50 


150 


6 


0.0675 


65 


4950 


14100 


37.5 


4.5 


50 


150 


7 


0.074 


60 


4100 


11350 


32 


5.5 


75 


200 



Table II 

Daily Generation of 10,26,34 Unit System 



Hon 


hi 


lOOait 


26 IB 


3* IB 




103 


1820 


1025 




100 


1500 


1000 




SO 


1720 


900 




85 


1700 


850 




103 


1750 


1025 




140 


1)10 


1400 




1S7 


2050 


1910 




240 


2400 


2400 




285 


2600 


2850 


10 


315 


2600 


3150 




330 


2620 


3300 


12 


340 


2580 


3400 


13 


32! 


2590 


3275 


14 


2S5 


2570 


2950 


15 


270 


2500 


2700 


16 


255 


2350 


2550 


17 


212 


2390 


2725 


1! 


320 


2480 


3200 


IS 


330 


2580 


3300 


20 


2S0 


2620 


2900 


21 


213 


2600 


2125 


22 


1(5 


2480 


1650 


23 


130 


2150 


1300 


24 


115 


1S0O 


1150 



The simulated demand set, corresponding generation 
schedule, the minimum operating cost and the 
computational time for the utility system is given in the 
Table III. The status of unit i at time t and the start-up / 
shut - down status obtained are the necessary solutions and 
are obtained for DP, LR, PSO, IPSO methods for utility 
system. 

Table III 

Optimal Generation Schedule For Utility System Satisfying 24 Hour 

Demand Along With Its Total Operating Cost 



It 



9 ii ii i: u it r h 17 is is :i :i ;: n ;i 



« ai i'i til tit a in tti a tti t'l tit m a 'it w '« +■• a tti t'l a tti u 



CI' <H <!>' II ill 11 ill ill tC V II' 1117' 



itt Ti ii ls.:i:i la'ii n n n n nun n n 
■ i i i i i i i i i in umiuim: 
171 i" is: si 191 in i«i w a w m n us is m 



i'i n 
U7: us: 1 

1H li] 



i n n ii? us 



a 1SI 13 1?) 191 
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performance of PSO 




20 25 30 
Number of Iteration 



40 45 



50 



Figure 3. The normalized cost for an optimal generation schedule versus 
the number of iterations of the PSO operation 

The comparison of the total costs and Central Processing 
unit (CPU) time is shown in Table IV for utility system, 10, 
26, and 34 generating unit power systems. The demand for 
24 hour time horizon is just simulated and it is not the 
actual demand which is practically satisfying by the units. 
An optimal generation schedule for each period and the 
total operating cost for the whole 24 periods are obtained. 
The IPSO algorithm contributes in determining the optimal 
generation schedule for a particular demand. The 
performance of PSO for a particular demand is depicted in 
the "Fig. 3". 

Table IV 
Comparisons of cost and CPU time for Utility & IEEE systems 



S^tem 


Methods 


Total Cost 
(p.u.) 


C PU Tim? 
(Sec) 




DP 


1 00000 


605 


7 - Unit Utility System 


LR 
PSO 


96481 
92690 


578 
519 




IPSO 


0.92420 


510 




DP 


1 00000 


325 


1 -Unit EEE System 


LR 
PSO 


0.94123 
90905 


279 

218 




IPSO 


0.90710 


211 




DP 


1 00000 


509 


26 - Unit EEE System 


LR 
PSO 


0.95968 
0.91273 


495 
440 




IPSO 


0.90S02 


431 




DP 


1 00000 


1452 


34 -Unit EEE System 


LR 
PSO 


0.99910 
0.95790 


1368 
1320 




IPSO 


0.95185 


1307 



Given a demand, the PSO generates an optimal unit 
commitment with minimum cost. In Figure 3, the 
improvement of PSO is illustrated in terms of offering the 
commitment of units with minimum cost. The affixed 
graph is obtained for solving the power demand of 400 
MW by the seven unit's utility system. In every number of 
iteration, the cost of the schedule offering by the IPSO gets 
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reduced. After a certain number of iterations, the cost 
remains constant for all the remaining iterations, which 
means that there no more generation schedule is available 
with cost which lesser than the previous cost. For the 
evaluation of performance, we have solved the UCP by 
IPSO only and thus we have compared the total production 
cost and computational time taken by the proposed 
approach and by the PSO to solve the problem and found 
that it was reduced for the IPSO method for all the systems. 

Conclusion 

The proposed approach IPSO has performed well in 
solving the UCP by recognizing the optimal generation 
schedule. The approach has been tested for the seven unit's 
utility system with the consideration of load balance and 
spinning reserve constraints, which are the most significant 
constraints. Prior to test the system, we have trained the 
network by different possible combinations of the demand 
set and its corresponding optimal schedule using the BP 
algorithm. 

For the test demand set which consists of demand for 24 
periods, the hybrid approach effectively yields optimal 
generation schedule for the periods. In comparison with the 
results produced by the referenced techniques (DP, LR, 
PSO), the IPSO method obviously displays a satisfactory 
performance. There is no obvious limitation on the size of 
the problem that must be addressed, for its data structure is 
such that the search space is reduced to a minimum; No 
relaxation of constraints is required; instead, populations of 
feasible solutions are produced at each generation and 
throughout the process. 
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